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ABSTRACT 

Models of Galactic 1.275 MeV emission produced by the decay of the radionuclide 
22 Na have been computed. Several frequency-spatial distributions of novae have been 
investigated using recent results of nova rates and spatial distributions of novae in our 
Galaxy. These models allow us to estimate the lower-limit of the 22 Na mass ejected per 
ONe nova detectable with the future spectrometer (SPI) of the INTEGRAL observa- 
tory as a function of the frequency-spatial distribution of ONe novae in the Galaxy. 
Calculations using recent estimations of the expected 22 Na mass ejected per ONe nova 
show that the detection of the Galactic emission of 1.275 MeV photons will be difficult 
with the future spectrometer of the INTEGRAL observatory, whereas the cumulative 
emission around the Galactic center has some chances to be detected during the deep 
survey of the central radian of the Galaxy. 

Key words: gamma-rays: observations - novae, cataclysmic variables - white dwarfs 
- Galaxy: structure 



1 INTRODUCTION 

Observations with several instruments have provided upper- 
limits to the 1.275 MeV flux, related to 22 Na decay from 
t he G alaxy or from individual novae, (see Clayton & Hoyle 



(1974) for the original idea to look for Na lines in novae) 



Observations with the HEA03 7-ray spectrometer provided 
an upper limit of 4 x 10 -4 photons cm -2 s _1 on the Galac- 
tic accumulated flux at 1.275 MeV (Higdon & Fowler (p.987p, 



hereafter HF87). With the 7-ray s pectr ometer of the Solar 
Maximum Mission, Leising et al. (11988 ) found a 99% con- 
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fidence limit of 1.2 x 10~ photons cm - s~ on a steady 
1.275 MeV fl ux from the Galactic center direction. Iyudin 
et al. (1995), using observations of single novae, with the 



COMPTEL instrument on board the Compton Gamma- 
Ray Observatory, estimated a 2a upper-limit of 3 x 1CF 5 



photons cm s for any neon-type nova in the Galactic 



disk, which has been translated into an upper-limit of the 
ejected 22 Na mass of 3.7 10" 8 M . 

The next generation of 7-ray spectrometers will have 
the sensitivity required for the detection of the 1.275 MeV 
emission from classical ONe novae. In particular, SPI, the 
future spectrometer of the INTEGRAL observatory will be 
operational in 2001. With its high energy resolution (AE/E 
wO.2%), the spectrometer is designed for the detection of as- 
trophysical 7-ray lines in the 20 keV - 8 MeV energy range. 



The detection plane is made of 19 germanium (Ge) detectors 
and is surrounded by an active shield made with bismuth 
germanate (BGO) scintillators, in order to reduce the in- 
strumental background. SPI will be able to perform images 
with an angular resolution of «3° by using a coded aperture 
system made of tungsten elements. The mean fully coded 
field-of-view is «16° (for a more detai led description of the 
s pectr ometer see Mandrou et al. (1997) and Teegarden et al. 
( 1997 )) . Using c alcula tions of 22 Na yields in ONe novae by 
Jose & Hernanz ( 1995 ) and instrumental background predic- 
tions, Hernanz et al. ( 1997 ) and Gorhez-Gomar et al. (199c) 
estimated that the 1.275 MeV line from an individual ONe 
nova could be detected by SPI if its distance is less than « 
0.5 kpc (standard CO novae produce 22 Na in much smaller 
amounts than ONe novae). With a new analys is of the nu- 
clear reaction rates involved in 22 Na synthesis ( Jose, Coc & 
Hernanz 1999) a bit large r distance is obtained (around 1.5 
kpc, Hernanz et al. (1999)). 

The total Galactic flux at 1.275 MeV depends not only 
on the amount of 22 Na ejected per outburst but also on the 
distribution and the rate of Galactic ONe novae. However, 
the Galactic nova rate is poorly known, independently of 
the underlying white dwarf composition, because the inter- 
stellar extinction prevents us from directly observing more 
than a small fraction of the novae that explode each year. 
Several methods have been used to estimate the nova occur- 
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rence rate. Estimations based on extrapolations of Galactic 
nova observatio ns sug gest a rate in the range 50-100 yr" 1 (see 
Liller & Mayer (1987) and references therein). Observations 
of novae in other galaxies have revealed a correlation be- 
tween the nova rate and the infrared luminosity of the parent 
galaxy. However, the occurrence rates of Galactic novae com- 
puted by scaling the nova rate measured in external galaxies 
are low er than «50 yr -1 . With this method, Ciardullo et al. 
(1990) predicted a n ova r ate in the range 11-45 yr -1 , and 
Delia Valle & Livio ( |l994| ) a rate ^20 yr" 1 . This is signifi- 
cantly lower than estimatio ns bas ed on Galactic nova obser- 
vations. Recently, Shatter (1997) reconciled this difference 
by recomputing the nova rate with the Galactic nova data. 
He extrapolated the global nova rate from the observed one, 
accounting for surface brightnesses of the bulge and the disk 
components and correction factors taking care of any obser- 
vational incompleteness. With this method, Shatter (1997 
estima ted the nova rate to be 35±11 yr" 
( [1997| ) found a similar value (41±20 yr" 1 
Carlo technique with a simple model for the distribution 
of dust and classical novae in the Galaxy. As only ONe no- 
vae are important producers of 22 Na, the Galactic ONe nova 
rate is the relevant one for this work. The Galactic ONe nova 
rate is obtained by multiplying the total nova rate by the 
fraction of novae that results from ther monu clear runaways 
on ONe white dwarfs. Williams et al. dl985| ) deduced from 



. Hatano et al. 
) using a Monte- 



observations of abundances in nova ejecta a proportion of 
ONe novae from 20 to 40 percent. With a detailed model of 
white dwarf mass distribution in cataclysmic binaries, Ritter 
et al. (1991) estimated that 25-57 percent of observed nova 
outbursts should occur on ONe white dwarfs. Livio & Tru- 



ran (1994) reestimated the frequency of occurrence of ONe 



novae, in light of observations of abundances in nova ejecta. 
They concluded that, of the 18 classical novae for which de- 
tailed abundance analyses were available, only two or three 
had a large amount of neon and were ONe novae, whereas 
three other novae showed a modest enrichment in neon, cast- 
ing doubt on the type of the underlying white dwarf. Under 
these considerations, they estimated a fraction of Galactic 
ONe novae between 11 and 33 percent. 

In addition to the 22 Na mass ejected per outburst and 
the ONe novae rate, the distribution of ONe novae in the 
Galaxy is an important parameter that tunes the total 
1.275 MeV flux arriving on Earth. Some 22 Na mass upper- 
limits, resulting from 1.275 MeV total Galactic flux upper- 
limit measurements, are underestimated since they have 
been computed assuming that all the ONe novae are em- 
bedded in the Galactic Center. Indeed, a larger 22 Na mass 
can be ejected by ONe novae leading to a 1.275 MeV to- 
tal Galactic flux below the measured upper-limit if the ONe 
nova distribution is more diluted in the Galaxy. 

In this paper, we model the Galactic cumulative emis- 
sion at 1.275 MeV as a function of the ONe novae rate and 
the mean 22 Na yield per outburst, for several spatial distri- 
butions of novae in the Galaxy. The lower-limit for the mean 
mass yield of 22 Na ejected per nova, in order to detect the 
1.275 MeV cumulative emission, is predicted for the future 
observations of SPI, using the estimation of the sensi tivity 
of the futur e spec trometer developed by Jean ( 199£ ) and 
Jean et al. (1997). The method for the simulation of the 
Galactic 1.275 MeV emission and its results are described in 
section 2, whereas the method of observation with SPI and 



its application to the previous simulations are presented in 
section 3. Discussion and conclusions follow (section 4). 



2 SIMULATION OF THE GALACTIC 
1.275 MeV EMISSION 

2.1 Method 

The 1.275 MeV emission from Galactic ONe novae is cal- 
culated with a Monte-Carlo simulation. A simulation of a 
set of ONe novae that is representative of what could be 
the Galactic novae distribution at a given time, i.e. a given 
novae frequency-spatial distribution, is made (this set of no- 
vae is further called a 'Galaxy-test'). The two nova popula- 
tions (bulge and disk) have been considered. The population 
origin (bulge or disk), the position in galacto-centric coor- 
dinates (depending on the chosen population) and the age 
of ONe novae are chosen randomly according to the appro- 
priate distributions. The local-galactic coordinates (l,b) of 
each nova as well as its distance to the Sun are calculated 
assuming a Sun to Galactic center (hereafter GC) distance 
of 8 kpc. The 1.275 MeV flux (F 22 in photons cm" 2 s" 1 ) has 
been estimated using the relationship 



F 22 = 3.1710 



_ t 

5 M 22 e T 22 

A 22 T 22 d? 



(1) 



where M 22 is the 22 Na mass ejected per nova (in solar 
masses), t the age of the nova (in years), d its distance to 
the Sun (in kpc). A 22 and T22 are the mass number and 
the lifetime of 22 Na (3.75 yr), respectively. Limitation of the 
flux due to opacity of the ejecta has not been taken into 
account, si nce it is significant only in the first week after the 
explosion (Gomez-Gomar et al. 1998), which is a short time 
as compared to the 22 Na lifetime. The number of simulated 
ONe novae depends on the global (CO and ONe) nova fre- 
quency (/n), the total period during which ejected 22 Na is 
an effective emitter (N T T22), the proportion of novae in the 
bulge (pi,) and the proportion of ONe novae in the bulge 
(with respect to all novae in the bulge, qb) and in the disk 
(with respect to all novae in the disk, qd). The number of 
simulated ONe novae is calculated as 



N„ = /„ N T t 22 [p b q b + (1 - p b ) q d ] 



(2) 



It has been assumed that the emissivity of novae older than 
5t22 has a negligible contribution to the diffuse Galactic 
1.275 MeV emission; therefore, N T =5. The proportion of 
ONe novae in the bulge (with respect to all ONe novae) 
is given by 



bulge 
PONe 



Pb qb 



Pb qb 



(1 



Pb) qd 



(3) 



The simulation of a Galaxy-test provides a set of 
1.275 MeV line fluxes Fi for each nova i localized at U, hi. 
Distribution of the intensity and of the number of novae has 
been calculated with 3° by 3° size pixels, which correspond 
to the angular resolution of SPI. 



2.2 Models 
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2.2.1 Distribution of novae in the Galaxy 

For the purpose of this work, a Galaxy model with a disk 
and a spheroid is convenient to simulate the distribution of 
Galactic nova events. Several laws for the spatial distribu- 
tion in the disk and in the 'spheroid' (representative of the 
bulge) have been proposed in the literature. We have se- 
lected four models that differ significantly from each other. 
The fi rst and older of them is described in Higdon & Fowler 
The radial scalelength of the disk is derived from 



the sta rlight su face brightness distribution (de Vaucouleur; 
& Pence 1978) and the height dependence is from Pat 



terson ( |1984| ). HF87 use the analytic dependence for the 
spheroi d ligh t distribution from Bahcall et al. (1982) and 
Young ( 1976 ), also called the 'R 1 / 4 ' law (R is the distan ce to 
the GC), originally formulated by de Vaucouleurs ( 1948| ). By 
scaling the nova rate in the bulge of M31 to the bulge of our 
Galaxy, HF87 estimate a proportion of novae in the spheroid 
of 0.348. The second model has been used by Hatano et al. 



(1997) to estimate the spatial distribution and the occur- 
rence rate of Galactic classical novae. It is based on a model 



of the distribution of type la supernovae (Dawson & John- 
son 199^). The proportion of novae that occur in the bulge 
is set to 0.111 on the basis of an estimate of the bulge to 
total galaxy mass ratio. The third model is derived from the 
Galactic survey of the Spacelab InfraRed Telescope (IRT) 
that provides a reliable tracer of the distri butio n of G and 
K gia nt stars (see Kent, Dame & Fazio (1991) and Ken t 



(1992)). It has also been used by Prantzos & Diehl (1996) 



to estimate the contribution of old population stars (novae 
and AGB) to the Galactic 26 A1 emission at 1.8 MeV mea- 
sured by COMPTEL. Using the total infrared luminosity of 
the bulge and the disk, the derived proportion of novae oc- 
curring in the bulge is 0.179. The last model is ta ken fr om 
Van der Kruit 



199q ). It has been used by Shatter jl997| ) to 



estimate the nova rate in our Galaxy. This author assumes 
that the nova distribution follows the brightness profile of 
our Galaxy. Under this assumption, the proportion of bulge 
novae is 0.105. 

The adopted distributions are shown in Figures |lj ^ 
and ^[ which represent the probability to find a nova in the 
Galactic bulge at spherical radius R in a volume element 
(kpc - ' 5 , Figure Q), and the probabilities to find a nova in 
the Galactic disk at cylindrical radius p in an area element 
(kpc" 2 , Fi gure B) and at height z above the Galactic plane 
in a line element (kpc -1 , Figure ^). Table [l] summarizes the 
characteristics of the adopted models. 



2.2.2 Galactic ONe nova rate 

As presented in section 1, methods used to estimate the oc- 
currence rate of novae are based on extragalactic and Galac- 
tic nova observations. The first method can underestimate 
the Galactic nova rate since the nova count can be biased by 
the extinction of the Galactic disk. On the other hand, es- 
timations obtained with Galactic nova observations depend 
on the Galaxy model used for the extrapolation. The l ast es - 
timations of classical novae rates are t hose by Shatter ( 1997 ) 
(35±11 yr -1 ) and Hatano et al. (|l997| ) (40±20 yr -1 ). There- 
fore, for the calculations presented in section 2.3, we have 
adopted frequencies of Galactic novae ranging from 20 yr -1 
to 60 yr -1 and a proportion of ONe novae from 10 % to 



Figure 1. Radial distributions of novae in the Galactic bulge for 
four models (in kpc -3 ). HF87: model used by Higdon and Fowler 
(1987) derived from the spheroid light distribution. This last 
model has also been adopted by Van der Kruit, 1990 (VdK90). 
KDF91: model from Kent, Dame and Fazio (1994) correspond- 
ing to the distribution of G and K giant stars provided by the 
InfraRed Telescope observations. DJ94: model from Dawson and 
Johnson (1994) derived from the distribution of type la super- 
novae. This last model has been truncated at 3 kpc (see table M 
for details). 




Figure 2. Radial distributions (cylindrical coordinate) of novae 
in the Galactic disk (in kpc -2 ) for the four models discussed in 
figure [j] (see text for details). 



Figure 3. Height above the Galactic plane distributions (cylin- 
drical coordinate) of novae in the Galactic disk (in kpc -1 ) for the 
four models discussed in figure [l] (see text for details). 
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Table 1. Nova spatial distributions used for the simulation of the Galactic 1.275 MeV emis- 
sion. R is the distance to the GC, z is the distance perpendicular to the Galactic plane and p is 
the galactocentric planar distance. The distance from the GC to the Sun is Rq=8 kpc, n s and 
rid are the normalization factors for the spheroid and the disk respectively (in kpc -3 ). The 
proportions of novae in the spheroid, p b , are 0.348 (HF87), 0.105 (VdK90), 0.179 (KDF91) 
and 0.111 (DJ94) 

Model HF87: p h =3.5 kpc and 2 o =0.106 kpc. 
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Model VdK90: p h =5.0 kpc and z h =0.30 kpc. 
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Model KDF91: p^=3.0 kpc and 2^=0.170 kpc. Ko is the modified Bessel function. 
Disk n(z,p) = n^e *h <>h 

Bulge n(R) = n s 1.04 x 10 6 (g-f^j ) ~ 1 ' 85 i? < 0.938fcpe 
= n s 3.53 K (0^7) B> 0.938fcpe 

= > 5fcpc 

Model DJ94: Ph=5.0 kpc and 2 h =0.35 kpc. 

Disk n(z, p) 
Bulge n(R) 



n d e 'h 

H 3 +0.343 




R < 3fcpc 
i?, > 3kpc 



Tru- 



30 %, in agreement with the recent estimate (Livio 
ran 19i 4|). These values correspond to a lower and upper 



limit of the ONe nova rate of 2 yr 1 and 18 yr 1 
tively. 



respcc- 



large as 10 8 Mq for a 1.25 Mq ONe nova in the most fa- 



vorab le case of nuclear reaction rates ( Jose, Coc & Hernanz 



199E) 



2.2.3 22 N a yield range 

The most recent upper-limit for the 22 Na mass ejected by 
ONe novae is provided by Iyudin et al. (1995), who de- 
rived 3 10 _8 Mq for any neon-type nova in the Galactic disk, 
at a 2 tr level with COMPT EL. Recent theoretical estima- 
tions (Jose & Hernanz 1998) provide typical 22 Na yields of 
r;2 10 a M per ONe nova, well below such observational 
upper limit. However, the average yield of 22 Na ejected per 
nova is not yet well known. It is sensitive to the composition 
of the underlying white dwarf and to the degree of mixing 
between the core and the accreted enveloppe (which con- 
strain not only the ejected envelope mass but also the max- 
imum temperature reached during the runaway) a nd also 
to the convection efficiency (Starrfield et al. 1996), (Jose 



&i Hernanz 1998). It is important to remind that in some 
cases theoretical models predict total ejected masses which 
are smaller than some observed ones (which in turn have of- 
ten large determination uncertainties). Moreover, some un- 
certainties of nuclear reaction cross-sections affect the com- 
puted 22 Na yields: a net increase in the 22 Na ejected mass by 
a typical ONe nova is obtained, with a very recent update of 
the nuclear reaction network affecting 22 Na synthesis. The 
new 22 Na yields for typical ONe novae are larger by factors 
ranging between ~2 and ~10, resulting in ejected masses as 



2.3 Results: maps of the 1.275 MeV emission 

Examples of Galaxy-tests for the distributions of novae in 
the Galaxy presented in section 2.2.1 are shown in Figures ^, 
|E| ^|and For these simulations, the adopted Galactic nova 
rate is 35 yr -1 (Shafter 1997) and the proportion of ONe 
novae is 30%. As already shown by HF87, the 1.275 MeV 
intensity is irregularly distributed in the Galaxy. We have 
recovered the 22 Na upper-limit mass found by HF87 with 
the se nsitivity of HEAQ3 - 4 10 -4 photons cm -2 s _1 at a 3a 



level, (Mahoney et al. 1982) - by simulating their frequency- 



spatial distribution of Galactic novae (rate 46 yr _1 and spa- 
tial distribution as shown in Table |l|) and by integrating the 
flux between 330° and 30° longitude. 

The statistical distributions of the 1.275 MeV fluxes 
from the 12.5° around the GC and from the brightest nova 
of the Galaxy have been calculated using a large number of 
Galaxy-tests. The results are presented in figure ^ and ^ for 
an ONe nova rate of ~10 yr -1 and for the four models of 
spatial distribution. The KDF91 model provides the high- 
est mean flux in the GC region, whereas the distribution of 
fluxes emitted by the brightest nova are similar in the four 
models, since the distributions of novae distances to the Sun 
are similar up to ~7 kpc for all the models. 

Table | shows the mean 1.275 MeV flux from the GC en- 
tering in the SPI field-of-view normalized to the 22 Na mass 
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Figure 4. Example of modeled diffuse 1.275 MeV emission gen- 
erated by a Monte Carlo simulation for the HF87 model. The 
fluxes have been normalized with the 22 Na ejected mass (photons 
s _1 cm -2 Mq" 1 ). The relative intensity is plotted with 3° by 
3° pixels in a map in galactic coordinates. For this example of 
Galaxy-test, the adopted ONe nova rate is f^lO yr _1 . 



Model: VdK90 



Figure 5. Same as figure 111 but for the VdK90 model. 



ejected per year in the Galaxy. The uncertainties give an 
idea of the intensity fluctuation magnitude due to the lim- 
ited number of effective emitting novae in the SPI field-of- 
view. The normalized fluxes are directly related to the spa- 
tial characteristics of the chosen Galactic ONe novae distri- 
bution (see Table [j]). The total 1.275 MeV normalized flux 
of the Galaxy is also presented in Table |^. The intensity in 



Model: KDF91 



Figure 7. Same as figure Q but for the DJ94 model. 



Table 2. Mean simulated 1.275 MeV flux emitted by ONe novae 
in the 12.5° around the GC region and in the whole Galaxy. The 
flux Fjv is normalized to the 22 Na mass ejected per outburst and 
the rate of Galactic ONe novae (photons cm -2 s _1 Mg _1 yr). 
Values correspond to the four adopted models. 



Models F N (GC) F N (Total) 
HF87 
VdK90 
KDF91 
DJ94 



78.0 ± 2.4 507 ± 85 

56.8 ± 3.5 580 ± 287 

105.1 ± 2.9 508 ± 62 

56.4 ± 1.4 496 ± 119 



the GC region represents 15%, 10%, 21% and 11% of the 
total flux in the Galaxy for the HF87, VdK90, KDF91 and 
DJ94 models, respectively. 

We have investigated the variation of the normalized 
1.275 MeV intensity in the 12.5° around the GC as a function 
of the scaleheight and scalelength of the novae distributions 
in the case of a zero proportion of novae in the bulge (see 
below). For a 0.3 kpc scaleheight, the normalized intensity 
decreases from «89 to «31 photons cm~ 2 s _1 M0 _1 yr when 
the scalelength increases from 2.5 to 6.0 kpc. Although the 
novae distances to the Sun are decreasing with increasing 
scalelength, the distribution around the GC is more diluted 
and the flux entering in the field-of-view is therefore lower. 
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Figure 6. Same as figure 111 but for the KDF91 model. 



Figure 8. Histogram of the normalized flux (flux per 22 Na ejected 
mass per outburst) in the 12.5° around the GC obtained from 
1000 simulated Galaxy-tests. The four models of spatial distribu- 
tion are presented. 
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Figure 9. Histogram of the normalized flux (flux per 22 Na ejected 
mass per outburst) emitted by the brightest ONe novae of the 
Galaxy. The four models of spatial distribution are presented. 
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Figure 10. Example of the total 1.275 MeV normalized flux in 
the 12.5° around the GC as a function of time. This results from 
a simulation with the HF87 model and an ONe nova rate of 11.5 
yr -1 . The flux in the GC can be dominated by a single close nova. 



3 ANALYSIS OF THE GALACTIC 1.275 MeV 
EMISSION BY SPI 

3.1 Types of detection 

Each one of the above mentioned Galaxy-tests could be ob- 
served by a 7-ray spectrometer at a given date. This is true 
if the observation time is short as compared with the occur- 
rence period of ONe novae and with the time variation of the 
22 Na emissivity. These conditions are satisfied since the com- 
mon observation time of the future 7-ray spectrometer SPI 
will be ~11 days, that is shorter than both the average ONe 
novae occurrence period («45 days) in the whole Galaxy, 
and the decay period of 22 Na (3.75 yr, implying 1.275 MeV 
flux reduction of only 0.8% in 11 days). However, for an ob- 
servation time longer than few times the ONe novae occur- 
rence period, the time variation of the 1.275 MeV Galactic 
emissivity has to be accounted for. Figure [lo] shows an ex- 
ample of the total 1.275 MeV flux variations in the 12.5° 
around the GC. It results from a Monte Carlo simulation 
as described in the previous section. The flux has been nor- 
malized to the 22 Na ejected mass per outburst. In figure |Io| 
the flux in the GC is dominated by a single nova (which 
contributes to 91% of the total flux) at 8.6 yr. After 20 yr, 
the flux shows fluctuations with a maximum amplitude of 
±40% of the mean flux. Figure [H] shows the statistical dis- 
tributions of the relative variations between two epochs of 
the total 1.275 MeV flux in the 12.5° around the GC. They 
have been calculated using the flux variation of figure llOl 
Several delays between the two epochs are presented. For 
two observations 6 months apart, there are probabilities of 
15% and 21% to get a flux variation of more than 10% and 
less than -10% respectively. Even if the flux variations are 
not large, the weak additional spots that could appear in 
the field-of-view will modify the distribution of the pixel- 
intensity. This could modify significantly the analysis in the 
case of the detection of excesses in the spatial distribution 
(see section 3.1.2). 

As stated above (section 2.1), the number of novae and 
the distribution of the intensity have been calculated with 3° 
by 3° pixels. Therefore, several novae can be located in the 
same pixel (this happens mostly in the direction of the GC 
region) . However, the flux in those pixels can be dominated 
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-0.5 0.0 0.5 

relative flux variation (5f/f) 

Figure 11. Distributions of the relative variations of the 
1.275 MeV normalized flux in the 12.5° around the GC between 
two observations separated by delays of 0.02 yr, 0.5 yr, 2 yr, 7.7 yr 
and 19 yr. These distributions are obtained with the light curve 
presented in figure [l(j. 



by only one nova. In this case, the measured intensity can 
not be considered as a cumulative effect. Therefore, we state 
for a seldom-nova criterium when a single nova contributes 
to more th an 95 % of the total flu x in a pixel. 

Jean (1996) and Jean et al. (1997) have estimated the 
narrow line sensitivity of SPI for an on-axis point source 
by computing instrumental background and detection ef- 
ficiency with Monte-Carlo simulations includ ing detailed 
physics. However, Gomez-Gomar et al. (1998) have shown 
that the 1.275 MeV line should have a typical width of 20 keV 
(FWHM). With these conditions, the sensitivity is 2.2 10" 5 
photons cm -2 s _1 



at a 3(7 level and for 10 s («11 days) 
of observation time, instead of 6 10 -6 photons cm -2 s . In 
the present case, the Galactic ONe novae emission is dif- 
fuse and the spectrometer will also provide a distribution of 
intensity in pixels. 

Therefore, three types of detection have been investi- 
gated for the analysis of any Galaxy-test. 
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3.1.1 Detection of the cumulative emission at 1.275 MeV 
( case 1 ) 



The number of sigma for each pixel is calculated by com- 



The to tal 1.275 MeV flux impinging on the SPI dctcctora 

will generate counts that will be added to the instrumental 
background ones. It will be possible to estimate these source 
counts by substracting the instrumental background compo- 
nent from the overall detector rates without using imaging 
deconvolution. Contrary to the imaging method, where the 
background rate is estimated using detectors that are not 
illuminated by a source (i.e. half of the total number of de- 
tectors - which are below opaque mask elements), the instru- 
mental background is estimated using either observations 
pointing out of the galactic plane - i.e. the on/off method 
used in balloon-borne spectrometers - or background mod- 
elling. The latter can consist of estimating the instrumen- 
tal background under the line of interest either by using 
data obtained during the observation (rate around the stud- 
ied line, other background feature rates), or by fitting the 
background temporal variations before and after the obser- 
vations, or both of the two methods. 

For this study, we simply consider the measurement of 
the cumulative emission with the on/off method. The sensi- 
tivity is estimated using the calculation of Jean (1996) (see 
also SPI - Science Performance Report ( ^PI - Science Per- 



forman ie Report 1996)) assuming that the time spent to 
measure instrumental background is equal to the GC ob- 
servation time and that there is no more modulation of the 
source photons by the coded aperture. In this case the expo- 
sure is on average equal to the point source exposure since 
half of the flux irradiates the detection plane (the other half 
being hidden by the coded mask elements). Since the counts 
recorded over all the detectors are taken into account, the 
background is increased by a factor «2 with respect to the 
imaging case. Therefore, since the sensitivity increases with 
the square root of the background, the sensitivity of such a 
mode of detection is roughly the sensitivity for an on-axis 
point source (in the imaging case) multiplied by the square- 
root of 2 (i.e. 3.1 10~ 5 photons cm -2 s _1 at a 3a level and 
for 10 6 s of observation time). In this method, the single 
information that will be available is the number of detected 
1.275 MeV photons coming from the region in the field-of- 
view of SPI. Although such a detection will give the proof of 
22 Na decay in the Galaxy, it will not provide any information 
on its spatial distribution in the observed region. However, 
galactic 22 Na yield can be estimated provided that assump- 
tions on the ONe distribution in the Galaxy are made. 



paring the flux with the sensitivity for a point source ( Jean 



et al. 1997). Afterwards, this number of sigma is converted 



in probabilities. The combination of the probabilities of the 
overall pixels allows us to estimate the probability that the 
observed distribution is not due to background fluctuation. 
Equation ^ shows an estimation of this probability P that 
is available when the significances in pixels are high enough. 
N is the total number of pixels in the SPI field-of-view, k 
is the number of pixels above a given significance threshold 
(2.5cr in this work) and Pi is the probability derived from 
the significance in the pixel i. 



P = 



(N -k)l 



(4) 



The distribution of the emission from the GC is consid- 
ered to be detected when the calculated probability leads to 
a significance that is above 3a for an observation time of 10 6 
seconds («11 days). A single nova close to the Earth can be 
the major source in the SPI GC field-of-view. Account for 
such an event makes a bias on the statistics of the cumu- 
lative 1.275 MeV flux. Therefore, to avoid such effect, the 
effective single nova pixel is removed from the analysis if its 
flux is more than 95% of the central 25° total flux. 



3.1.3 Detection of the brightest nova of the Galaxy in 
1.275 MeV (case 3) 

The flux from the brightest nova that can be observed from 
the Earth is also analyzed and considered to be detected if its 
value is above the sensitivity of the spectrometer SPI, for an 
observation time of ~11 days. The observation time needed 
to detect with SPI (at a 3a level) the brightest nova of the 
Galaxy-test has also been estimated. Equation (^) allows us 
to calculate the duration of the observation {Tza in days) 
to get a 3a detection of the flux / (photons cm -2 s _1 ) in a 



single pixel. 

T 3ct = 2.42 10" 5 f~ x + 5.28 10" 9 /" 2 



(5) 



We can recognize the well known variation of the sen- 
sitivity as a function of the square root of the observation 
time in the second term of the equation, that is dominant 
when the signal-to-noise ratio is low (i.e. low fluxes). The 
first term accounts for the statistical fluctuations induced 
by the counts from sources when the signal-to-noise ratio is 
not low anymore (i.e. for high fluxes). 



3.1.2 Detection of excesses in the spatial distribution of 
the intensity at 1.275 MeV (case 2) 

The total flux in the SPI field-of-view is distributed in sev- 
eral pixels. Some of the pixels can contain an intensity that 
leads to a significant signal. In this case, we obtain an ad- 
ditional information that is the location of the emission in 
the Galaxy. However, it is necessary to account for the re- 
duction of the significance of the detection with the number 
of pixels involved in the total flux. A rough estimation of 
the significance (number of a) of such a detection is derived 
by computing the probability that the observed intensity 
distribution in the field-of-view is due to background fluctu- 
ation. The lower the probability, the higher the significance. 



3.2 Probability to detect the 1.275 MeV line 
with SPI 

For a large number of Galaxy-tests, the probability to de- 
tect any type of emission can be estimated by calculating 
the fraction of time that the simulated 1.275 MeV flux is 
above the upper-limit of a given experiment (HEA03 or 
SPI). For the particular case of the future spectrometer of 
INTEGRAL, the average times of observation to obtain de- 
tections of the three types mentioned above have been de- 
rived, as well as the average number of novae that contribute 
to the diffuse emission of the GC region. Each Galaxy-test 
has been analyzed in order to check whether SPI would de- 
tect: the cumulative emission at 1.275 MeV in the 12.5° 
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around the GC region (case 1), an excess in the distribu- 
tion of the emission (case 2) and the brightest nova of the 
Galaxy (case 3), for a given average 22 Na yield per ONe 
nova outburst. The time needed for a 3a level detection is 
also computed. Two hundred Galaxy-tests have been simu- 
lated for each model and for a given 22 Na yield in order to 
calculate the number of Galaxy-tests that lead to a detec- 
tion by SPI with an observation time of 10 6 s. The number 
obtained, divided by the total number of Galaxy-tests, gives 
an estimation of the probability to detect the 1.275 MeV line 
from ONe novae with SPI. 

Figure [li] shows the mean observation times for the 
three types of detection with SPI, as a function of the av- 
erage 22 Na mass ejected per nova. The adopted ONe nova 
rate is 10 yr _1 and the spatial distributions of novae are those 
described in section 2.2.1. For low values of 22 Na yield, the 
observation times are very long and do not make sense since 
the limited duration of the INTEGRAL mission does not al- 
low for observation times larger than 2 years. Moreover, as 
the Galactic 1.275 MeV intensity and its distribution would 
change due to the decay of 22 Na and to the explosion of 
additional novae, the analysis of the intensity distribution 
(case 2) and of the brightest nova flux (case 3) are not valid 
anymore. However, other analyses (cases 1) can be possible 
with exposures in the limit of the allocated observation time. 

Figure [l^ shows the fraction of time that SPI would 
detect, in «11 days, an excess in the distribution of the 
emission from the GC region (case 2) as a function of av- 
erage 22 Na masses ejected per nova, for Galactic ONe nova 
rates ranging from 2 to 18 yr _1 (i.e., an estimation of the 
probability to detect this emission). In this analysis, it is 
not possible to normalize the 22 Na mass ejected by the ONe 
nova rate, because when the frequency of novae increases, 
the (l,b) distribution of intensity changes, resulting in an 
estimation of the significance that implies more pixels and 
not only more flux. Therefore, the sensitivity to an excess in 
the 1.275 MeV intensity distribution is not directly propor- 
tional to the amount of 22 Na ejected per year. 

Our results allow for the determination of the 22 Na 
yield lower-limit that can be detected by the spectrome- 
ter of INTEGRAL. We state that this lower-limit of 22 Na 
ejected mass corresponds to a probability of detection of 
90%. The lower-limit depends also on the nova rate, since 
more ONe novae lead to higher cumulative GC region flux, 
on the Galaxy model and on the type of detection. Taking 
the 22 Na mass values at 90% of probability in Figure [l^ we 
can calculate the lower-limit as a function of the Galactic 
ONe nova rate for case 2. Similar calculations have been 
performed for the detection of the cumulative emission in 
the SPI field-of-view (case 1) and of the brightest nova of 
the Galaxy (case 3). The results are displayed in Figures 
|l4| for the four adopted spatial distributions, as a function 
of the nova rate. As it has already been shown in Figure 
|l2] , it is clear that the future spectrometer SPI will be more 
sensitive to the brightest nova of the Galaxy and the total 
diffuse flux than to a detection of an excess in the distribu- 
tion of the emission. The 22 Na lower-limit masses needed to 
detect the brightest nova (case 3) are similar for the four 
models, since the distributions of novae distances to the 
Sun associated to these models are similar up to ~7 kpc 
from the Sun. Concerning the 22 Na lower-limit masses for 
the cumulative emission (case 1), it can be seen in Figure 



|l4| that models KDF91 and HF87 provide lower values than 
VdK90 and DJ94 models. The radial scalelengths for the 
former are 3.0 kpc and 3.5 kpc, respectively, whereas the 
models VdK90 and DJ94 have higher radial scalelengths, 
5 kpc in both cases (see Table ^). Since novae are closer to 
the GC in the KDF91 and HF87 models, their cumulative 
1.275 MeV flux is higher and more compact (see Figures ^j, 
|E| ^| and Q). For a given ONe nova rate, the 22 Na lower-limit 
mass decreases with decreasing scaleheight of the Galactic 
disk: the lower the z h (i.e., models HF87 and KDF91), the 
higher the cumulative 1.275 MeV flux at the GC (see Table 
^|) and the smaller the 22 Na lower-limit mass (see Figure |l^) . 
However, this variation is not significant: for instance, when 
the scaleheight increases from 0.1 kpc to 0.45 kpc the 22 Na 
lower-limit mass increases from ~7 10 -7 Mq to 9 10 -7 Mg, 
if we use the DJ94 model with variable scaleheight. More- 
over, the limited angular resolution of SPI constrains to a 
lower-limit the observable disk scaleheight. Indeed, in the 
case of an exponential disk model (without bulge contribu- 
tion), simulations show that the estimated 22 Na lower-limit 
mass does not change for zj, values lower than ~0.15 kpc and 
~0.25 kpc, which correspond roughly to a scale-angle of 3° 
for a Galactic disk radial scalelength of 5 kpc and 3 kpc, 
respectively. 

It is not clear whether the Galactic bulge should be a 
source of 1.275 MeV photons, since it has been proposed 
that it is 10 10 years old and white dwarfs responsible of 
ONe novae should be much younger (the zero age mass of 
the white dwarf progenitor being probably larger than ~10 
Mq). However, the ages of novae do not only depend on 
the main sequence lifetime of the progenitor but also on 
the time necessary for the onset of the Roche lobe overflow 
in the binary system and the time between the onset of 
accretion and the outburst phase. Taking into account the 
eventuality that there is a lack of ONe n ovae i n the bulge 
(see the recent paper by Delia & Livio (1998)), we have 
estimated the 22 Na mass lower-limits in the case of ONe 
novae only distributed in the Galactic disk. Figure |l5| shows 
the results of the simulations in this case for the four models 
adopted: the 22 Na mass lower-limits that SPI can detect are 
in general lower, since novae are concentrated in the disk 
and, therefore, they are closer to the Earth. 

As a final result, extracted from those presented in Fig- 
ure [l2| Figure [l(] shows the observation time needed to have 
90% of chance to detect with SPI the 1.275 MeV line emitted 
by the accumulated 22 Na from the GC region (case 1) as a 
function of the mean 22 Na mass ejected per ONe nova. These 
results are displayed for the four distributions of novae in the 
Galaxy. For the HF87 and KDF91 spatial distributions and 
our estimations of 22 Na yields, observation times compatible 
with those planned for the deep survey of the central Galaxy 
are obtained. Therefore, cumulative emission in the GC may 
be detected with SPI. 



4 DISCUSSION AND CONCLUSIONS 

It is important to remind that the total Galactic flux at 
1.275 MeV depends not only on the amount of 22 Na ejected 
per outburst and the ONe novae rate but also on their distri- 
bution in the Galaxy. The ejected 22 Na mass upper-limit de- 
rived by HF87 using HEA03 measurement was estimated for 
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Figure 12. Mean observation time needed to detect the cumulative (diffuse) emission in the GC (case 1), the distribution of the emission 
(case 2) and the brightest nova (case 3) at 1.275 MeV with the future spectrometer SPI, as a function of the mean 22 Na yield per nova. 
The error bars have a width of lc. The adopted Galactic ONe nova rate is 10 yr - x . Each point results from the processing of 200 
Galaxy- tests. Each panel corresponds to a different spatial distribution of novae (see text for details). 
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Figure 13. Fraction of time during which the distribution of the diffuse emission, in a sphere of 25° diameter around the GC would be 
detected by SPI (case 2), as a function of mean 22 Na yields per nova, calculated for ONe nova rates ranging from 2 to 18 yr _1 and for 
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Figure 16. SPI observation time necessary to have 90% of chance 
to detect the cumulative 1.275 MeV emission of the GC (case 1) 
as a function of the mean 22 Na mass ejected per nova and for 
several Galactic distributions. 
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Figure 17. Galactic ONe nova rate needed to have 90% of chance 
to detect the cumulative 1.275 MeV emission from the GC region 
(between 330° and 30°) with HEA03 as a function of the mean 
22 Na mass ejected per nova and for several Galactic distributions. 
The dot shows the estimation of HF87. 



only one distribution and an ONe nova rate of 11.5 yr -1 . Fig- 
ure [L7| shows the 22 Na mass upper-limits computed with the 
1.275 MeV upper-limit flux of HEA03 considering other pos- 
sible frequency-spatial distributions of novae in the Galaxy. 
It shows that the 22 Na mass upper-limit lies between 2 and 
10 icr 7 M per nova. 

The spectrometer of INTEGRAL will be more sensitive 
for the detection of the brightest nova of the Galaxy and 
of the cumulative diffuse emission of ONe novae in the GC 
region by an "on-off" type observation than for the deter- 
mination of their intensity distribution. But it is not clear 
whether the brightest novae will be first detected in the vis- 
ible due to the interstellar extinction. Nevertheless, long ex- 
posures in the Galactic plane region may allow for their de- 
tection with SPI. 

According to Jose, Coc & Hernanz (1999), the 22 Na av- 
erage yield could be between 3. 10 -9 and 1.2 10~ 8 Mo per 
nova. Therefore, there are few chances to detect the Galac- 
tic diffuse 1.275 MeV emission with SPI with only ~10 days 



of observation. However, 80 days of observation of the GC 
region could already give constraints to the mean 22 Na yield 
in ONe novae as a function of their frequency-spatial distri- 
bution in the Galaxy (see Figure |l6| ) . 

It is worth mentioning that the adopted width of the 



1.275 MeV line (20 keV) is quite large, corresponding to the 
fastest ONe novae, which have the largest ejection velocities 
(and also the largest ejected 22 Na masses). As the correla- 
tion between 22 Na ejected mass and ejection velocity is not 
well established (neither theoretically nor observationally) , 
and the braking of the ejected mass during the lifetime of 
22 Na is in principle small but not null, we are in the worst 
situation when we adopt such broad lines (corresponding to 
the largest velocities and without any braking effect). In the 
other extreme case of narrow lines (less realistic but interest- 
ing to make a comparison), the sensitivity of SPI is 6 10 -6 
photons cm~ 2 s" 1 at a 3a level and for 10 6 s (~11 days) 
of observation time, for a point source, which is ~3.7 times 



better than the one for our broad lines (see section 3T) and 
allows for the detection of lower fluxes by the same factor. 
Consequently, a positive detection of the 1.275 MeV line by 
SPI would be achieved either with a smaller 22 Na mass or 
with a smaller observation time. 

Gamma-ray observation of novae would provide infor- 
mation not only on their eruption mechanisms and the nu- 
cleosynthesis processes involved in their explosion but also 
on their distribution and their rate in the Galaxy (e.g. pro- 
portion in the bulge, scaleheight in the Galactic disk) since 
the problem of the interstellar extinction does not appear at 
this energy range. 
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